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ABSTRACT 

We present an analysis of archival Chandra observations of the mixed-morphology rem¬ 
nant 3C400.2. We analysed spectra of different parts of the remnant to observe if the plasma 
properties provide hints on the origin of the mixed-morphology class. These remnants often 
show overionization, which is a sign of rapid cooling of the thermal plasma, and super-solar 
abundances of elements which is a sign of ejecta emission. Our analysis shows that the thermal 
emission of 3C400.2 can be well explained by a two component non-equilibrium ionization 
model, of which one component is underionized, has a high temperature (kT sa 3.9 keV) 
and super-solar abundances, while the other component has a much lower temperature 
(kT ss 0.14 keV), solar abundances and shows signs of overionization. The temperature 
structure, abundance values and density contrast between the different model components 
suggest that the hot component comes from ejecta plasma, while the cooler component has an 
interstellar matter origin. This seems to be the first instance of an overionized plasma found in 
the outer regions of a supernova remnant, whereas the ejecta component of the inner region is 
underionized. In addition, the non-ionization equilibrium plasma component associated with 
the ejecta is confined to the central, brighter parts of the remnant, whereas the cooler compo¬ 
nent is present mostly in the outer regions. Therefore our data can most naturally be explained 
by an evolutionary scenario in which the outer parts of the remnant are cooling rapidly due to 
having swept up high density ISM, while the inner parts are very hot and cooling inefficiently 
due to low density of the plasma. This is also known as the relic X-ray scenario. 

Key words: ISM: supernova remnants - supernovae: general - supernovae: individual: 
3C400.2 


1 INTRODUCTION 

An important and not well-understood class of supernova remnants 
(SNRs) are the so-called thermal composite or mixed-morphology 
remnants (MMRs jRho & Petre| ( |1998| >; |Lazendic & Slane] ( |2006| l; 
|Vmk| ( |2012a) ). These remnants are characterised by thermal X-ray 
emission that is centrally peaked whereas the radio emission has the 
familiar shell-type morphology. They have several other interest¬ 
ing characteristics: they are often associated with GeV gamma-ray 
emission (e.g. | Giuliani et al.|201 0; Uchiyama et al. 2012}, there is 
evidence for enhanced metal abundances in some of them ( |Lazen-| 
|dic & Slane|2006| although they are usually mature remnants, and 
they often show spectral features associated with strong cooling 
in the form of radiative recombination continua (RRCs) or strong 
He— alLy—a X-ray line ratios of alpha-elements such as Si, S and 
Ar (Kawasaki et al.|2005} . 

The centrally peaked X-ray emission of MMRs poses a prob¬ 
lem for SNR evolution models. They typically have an age on the 
order of 20.000 years and are therefore expected to have a shell¬ 
like density structure with a hot, tenuous plasma in the centre, 


based on a Sedov evolutionary scenario. A flat interstellar matter 
(ISM) structure is therefore unsuited to create the centre-filled X- 
ray morphology that is observed, if the temperature is relatively 
uniform across the remnant. |Rho & Petre| ( j 1998} mention two pos¬ 
sible scenarios for the formation of the emission structure typical 
for MMRs, the relic X-ray emission and the evaporating cloudlet 
scenario i jWhite & Long||T99l} . The former scenario is the most 
simple scenario, in which the outer layers of the remnant have be¬ 
come radiative and have cooled strongly so that they hardly emit in 
the X-ray band, whilst the centre consists of hot plasma that has not 
yet cooled below 10° K. This scenario requires a high surrounding 
ISM density. [White & Long| l |l99l} give a self-similar solution for 
a supernova that exploded in an environment where it is surrounded 
by a number of small, dense clouds. Due to the large filling fraction 
and small size of the clouds they do not alter the dynamics of the 
forward shock, and they survive the passage of the forward shock 
due to their large density. The clouds then slowly evaporate due to 
heating by thermal conduction, which increases the density and de¬ 
creases the temperature in the centre of the remnant. Of these two 
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scenarios, the relic X-ray emission scenario was preferred by [Har-| 
rus et al.|(1997} to explain the morphology of W44. One year later, 
Cox et all | ( } 1 999| > ; | Shelton et al.| ( |1999| > used a similar scenario as 
Harms et aT](1997| > to model the characteristics of the MMR W44, 
but added a density gradient for the ISM and thermal conduction. 
Thermal conduction in the model by |Cox et al. 1(1999) smoothes the 
temperature gradient from the centre to the outer layers, thereby re¬ 
ducing the pressure in the centre. A lower central pressure reduces 
the need to expand, which allows for a higher density to remain 
in the centre of the remnant. A problem with the scenarios using 
thermal conduction is that it is not a priori clear whether it can 
be important in SNRs, since it is strongly suppressed by magnetic 
fields (Spitzer| 198 l)|Tao| 1995) . 

Clear evidence for overionization of thermal SNR plasmas is 
almost exclusively found in MMRs. Thermal plasmas in supernova 
remnants are often found in an underionized state, and reach ioniza¬ 
tion equilibrium on a density dependent timescale t ~ 10 12 ' 5 /n e 
s i jSmith & Hughes|2010) l. Due to their ages and their frequent as¬ 
sociation with high density regions, MMRs are expected to have 
a plasma that is in ionization equilibrium. The fact that there is 
ample evidence for overionization, means the cooling rate of the 
plasma was faster than the recombination rate. How the rapid cool¬ 
ing proceeds is still unclear. Efforts have been made to localise re¬ 
gions in MMRs which are cooling rapidly, as an association with 
a high density region may mean that thermal conduction is impor¬ 
tant. For example: |Miceli et al.|(2010) ; |Lopez et al.| ( |2013| > find that 
the amount of cooling increases away from a molecular cloud for 
W49B, suggesting adiabatic cooling as the dominant cooling mech¬ 
anism. |Uchida et al.| < (2012j > suggest the cooling timescale of the 
plasma cannot be explained by thermal conduction, and that adia¬ 
batic expansion is probably the dominant cause for cooling there. 
In addition, |Broersen et al.| ( |2011 | l have shown that even at small 
expansion rates, simple adiabatic cooling combined with cooling 
through X-ray radiation can lead to a cooling rate of the plasma 
larger than the recombination rate. So far observational evidence 
for strong thermal conduction is lacking. 

The ionization state of a thermal plasma can be determined 
by using both the electron temperature T e and the so-called ion¬ 
ization temperature T z (Masai |[T997) . The ionization temperature 
is the electron temperature that one would deduce from the ioniza¬ 
tion state of X-ray emitting elements alone, which may be different 
from the thermodynamic temperature in case the plasma is not in 
ionization equilibrium. When T z < T e , the plasma is underionized, 
when T z = T c the plasma is in collisional ionization equilibrium 
(CIE) and when T z > T e the plasma is overionized. There are dif¬ 
ferent ways in which T z can be determined. Kawasaki et al.|(2002) 
determine T z by comparing the observed Ly —a / He —a ratio of 
a certain element to the ratios of CIE plasmas of different temper¬ 
atures. The temperature of the CIE plasma that produces the ob¬ 
served ratio is then T z . A different method was used by|Ohnishi et| 
|ak](2011), who determined T z by using the CIE model in SPEX to 
fit the plasma. This model has an additional parameter £, compared 
to the XSPEC CIE model. This parameter can be used to simulate 
a non-equilibrium plasma so that T z = £ x T e . Finally there is 
the method used by |Broersen et al.|(20 1 1) > and |Uchida et al.| ( [2012| >, 
who fit a plasma using the SPEX NEI model with initial temperature 
kT\ > kTz. kT\ has a similar function as the ionization tempera¬ 
ture, but it has a slightly different physical meaning as it describes 
the initial ionization state. The NEI model then follows the ioniza¬ 
tion state as a function of ionization age as usual, assuming a rather 
sudden drop of electron temperature. Note that although using the 
He—a / Ly—a line ratios of alpha elements can show that a plasma 



Figure 1. The ionization temperature of a plasma versus the ionization age. 
The NEI plasma starts at an electron and ionization temperature of 4.0 keV, 
and it evolves after the electron temperature instantly drops to 0.2 keV. The 
figure illustrates that the NEI and line-ratio/ionization-temperature methods 
for characterizing overionized plasma do not provide identical ionization 
states. Moreover, the differences depend on the element: at certain ioniza¬ 
tion ages, different elements may have different ionization temperatures. 

is overionized, it has quite a large systematic error in determining 
the T z (Lopez et al.|2 013l. This is the result of the fact that there 
are emission lines of other elements which contaminate the He—a 
/Ly —a ratio. The strengths of the contaminating lines are T e depen¬ 
dent, so that the systematic error is T e dependent as well. This can 
be taken into account by calculating the CIE line ratios including 
the contaminating lines. 

The difference in ionization temperature between the line ratio 
and NEI model methods is displayed in Fig. [T] This figure shows 
the evolution of the ionization temperature as a function of ioniza¬ 
tion age, for an NEI plasma that is cooling abruptly from 4 keV 
to 0.2 keV. It is clear from the figure that different elements show 
different ionization temperatures as the plasma evolves, so that a 
plasma cannot be characterised by a single ionization temperature. 
This also explains why different elements in SNR plasmas often 
show different ionization temperatures, as found by |Lopez et al.| 
(e.g. |2013] l. On the other hand the NEI model is not a perfect model 
for a cooling plasma, as it assumes the electron temperature in¬ 
stantly drops to a certain value. In reality the electron temperature 
will drop more gradual, so that the plasma evolution will be slightly 
different as a result of different recombination and ionization rates. 

Here we present the first analysis of a Chandra X-ray observa¬ 
tion of 3C400.2 (also known as G53.6 - 2.2), which is an important 
member of the MMR class. This remnant has centrally peaked X- 
ray emission, as shown by Einstein IPC and ROSAT observations 
(Long et aL 1991; Sa ken et al.|l9 95 l. The radio morphology can 
be described as two overlapping circular shells of diameters 14’ 
and 22’ (Dubner et al.|199~4j >. This has led to the speculation that 
3C400.2 might be two supernova remnants in contact with each 
other, which would make it a rare event. [Yoshita et aL|(2001fr con¬ 
clude however, using ASCA observations, that the remnant is the 
result of a single supernova explosion, based on the similar plasma 
properties found in the two shells. Hydrodynamical simulations, in¬ 
cluding thermal conduction, show that the morphology of the rem¬ 
nant can be explained by a supernova exploding in a cavity, where 
the larger shell is a part of the remnant expanding into a lower den¬ 
sity region than the smaller shell (Schneiter et al.|[2006). This is 
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Figure 2. Rosat PSPC image in an inverted grey scale, with NRAO VLA 
Sky Survey contours overlaid, The Chandra ACIS-I field of view is indi¬ 
cated in red. 


consistent with HI observations performed by Giacani et al. [< [ 1998j> , 
who find a denser region to the northwest of the remnant where the 
smaller shell is located. In the optical, the remnant is characterised 
by a shell-like structure with a smaller radius (about 8’) than the 
radio shell ( [Winkler et al.|1993|lAmbrocio-Cruz et al.l2006| . The 
optical emission is located in regions of low X-ray emissivity. Op¬ 
tical emission suggest that those parts of the remnant contain ra¬ 
diative shocks, and are therefore cooling efficiently. Distance esti¬ 
mates to 3C400.2 range f rom 2.3 - 6.9 kpc (|Rosado|1983[ |Milne| 
|1979[|Dubner et al.| 1 994]|Giacani et al.| 1 998) . However, the dis¬ 

tance estimates based on radio data are obtained with the uncertain 
E — D relationship, and the kinematic estimate of |Rosado] l jl983[ ) is 
based on interferograms of a small part of the remnant. We there¬ 
fore consider the most recent distance estimate of 2.3 ± 0.8 kpc 
( |Giacani et al.| 1998] ), which is based on HI measurements, as the 
more reliable one and we will use a round value of 2.5 kpc for the 
distance throughout this paper. 

MMRs and in particular overionization of plasmas have gath¬ 
ered increasing attention over the past few years. In this work, we 
therefore aim to characterise the plasma properties of 3C400.2 in 
order to increase our understanding both of overionized plasmas, 
and the possible evolutionary scenarios of MMRs. We start with 
the data reduction and spectral analysis of different regions of the 
remnant, which is followed by a discussion on the measured plasma 
properties. We end with the conclusions. 


2 DATA ANALYSIS AND RESULTS 
2.1 Data reduction 

In this paper we report on the analysis of a 34 ks archival Chandra 
observation (ObsID 2807) taken on August 11, 2002 with Chan¬ 
dra observing in imaging mode with the ACIS-I CCD array. The 
total number of counts in the 0.3-7.0 keV band in the full spec¬ 
trum is ~ 58k. We extracted spectra using the default tasks in the 
Chandra analysis software CIAO version 4.5, using the task specex- 
tract to create spectra and weighted responses (RMF and ARF). 
We removed point sources from the spectral extraction regions. The 
spectra where grouped so that each bin contained a minimum of 15 


counts. The remnant covers the whole area of the ACIS-I chips (see 
[2j>, and there was no region available for background extraction. We 
therefore used the standard ACIS-I background files to create back¬ 
ground spectra. We used SPEX version 2.03 l |Kaastra et al.|1996| for 
the spectral modelling. The NEI model in SPEX is an extended ver¬ 
sion of the MEKAL model used in XSPEC ( |Amaud|1996| l, with the 
added advantage that the initial temperature of the NEI model can 
be varied to mimic an overionized plasma. 

2.2 Spectral analysis 

The Chandra ACIS-I field of view covers the part of the remnant 
that is mostly associated with the smaller radio shell. It includes 
the brightest part of the remnant in X-rays, as can be deduced from 
the ROSAT image in Fig. [2] With our spectral analysis, we aim to 
find differences in plasma properties between different parts of the 
remnant, and to investigate if this remnant conforms to the general 
properties of MMRs mentioned in the introduction, with regards to 
overionization and abundances. 

We fitted the spectra using absorbed non-equilibrium ioniza¬ 
tion (NEI) models, of which the parameters are the ionization age 
r = n e t, the electron temperature fcT 2 , the elemental abundances 
and the normalisation n e nnV. A CIE plasma is identical to an NEI 
plasma when r ^ 10 12 ' 5 cm -3 s. In addition to the above listed 
parameters, the SPEXNEI model has the initial temperature of the 
plasma, fcTi, as an optional parameter. As mentioned in the intro¬ 
duction putting fcTi > fcT 2 makes the NEI model mimic an ove¬ 
rionized plasma, where the ionization state of the plasma is deter¬ 
mined by fcTi and r, while the continuum shape is determined by 
the electron temperature fcT2 = fcT e . The method to check for ove¬ 
rionization using Ly — a and He — a line strengths is not feasible 
for 3C400.2, since elements with isolated emission lines such as Si 
and S show no significant Ly — a lines in this remnant. 

We aimed to first find a satisfactory fit for the region cover¬ 
ing almost the full area of the ACIS-I chips (see Fig. |4j, where 
we initially tried to fit the spectrum with an absorbed, single, un¬ 
derionized NEI model with fixed abundances, freeing abundances 
only when it led to a significant improvement in C-stat / degrees of 
freedom ( |Cash| 1979] >. The Ah was always allowed to vary. A sin¬ 
gle absorbed NEI model was not sufficient to obtain an acceptable 
C-stat / d.o.f. in any of the extracted spectra, however. 

The next step was to try to fit the spectrum with a double NEI 
model. Although the C-stat / d.o.f. improved significantly with re¬ 
spect to a fit with a single NEI model, the fit was still not acceptable. 
The only way the fit of the double NEI models became acceptable, 
was by allowing the fcTi of the cooler NEI component to vary, so 
that fcTi > fcT 2 , and therefore the NEI plasma is overionized. We 
show in detail why the overionized model fits the data better than 
an underionized model in section l2.2.1l 

Fig.[3]shows the spectrum of the full region. The best-fit model 
for this region consists of a fcT = 3.86 keV plasma with super¬ 
solar abundances plotted as a green dashed line, combined with a 
fcT 2 = 0.14 keV, overionized NEI component which is plotted as a 
blue dot-dashed line. The parameters of the model are listed in Tab. 
|T] It is clear from the figure that the high fcT model accounts for 
the bulk of the Fe-L (0.7-1.2 keV), Si XIII (1.85 keV) and S XV 
(2.46 keV) line emission, while the low fcT model accounts for 
the OVII-VIII (0.56 and 0.65 keV) line emission, and continuum 
emission in the form of O VIII, Ne X and Mg XI RRCs. The super¬ 
solar abundances of Si, S and Fe in the hot NEI component suggest 
an ejecta origin for the plasma. The abundances of the high fcT 
component as well as the emission measure have quite large er- 
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Figure 3. Spectrum of the region covering almost the whole area of the ACIS-I chips. The blue dash-dotted line shows the low kT overionized NEI component, 
while the green dashed line shows the high kT underionized NEI component. The parameters of the model are listed in Tab. [T] Indicated in the figure are 
the locations of the O VIII, Ne X and Mg XI RRCs which are present in the overionized model (blue dash-dotted line) and the locations of important emission 
lines. The green dashed line represents the hot NEI component. 



19:39:00 19:38:30 19:38:00 
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Figure 4. Chandra image in the 0.3 - 7.0 keV band with NVSS radio con¬ 
tours. The extraction regions of the spectra are labelled in yellow. 


rors. This is caused by the fact that the continuum shape at energies 
larger than 2.5 keV is not well defined by the data. An ill-defined 
continuum strength causes large formal errors in the abundances, 
since the strength of the continuum and the height of the abun¬ 
dances are anti-correlated. Nonetheless the abundances are signifi¬ 
cantly super-solar, which means that 3C400.2 belongs to the group 
of MMRs with super-solar abundances. The ill-defined continuum 


at higher energies may be the result of the uncertainty introduced by 
the blank-sky background subtraction. Overionization is found only 
in the lower kT plasma, and not in the high kT plasma. Making the 
initial temperature a free parameter in the hot NEI component did 
not improve the fit. There is therefore no evidence for overioniza¬ 
tion in the ejecta component, which is further corroborated by the 
absence of strong Ly — a lines of Si and S. The abundances of the 
low kT NEI component are sub-solar in the case of NE and O and 
solar for all other elements, which suggest a swept-up ISM origin 
for this plasma. 

Overall the fit is acceptable at a C-stat / d.o.f. of 257 / 227, 
although there are some significant residuals. The most notable 
residual feature is found at the position of the SixiII-XIV line at 
1.85 keV. The model fits this line partly with the Mg XI RRC from 
the cool NEI component, and partly with Si emission from the hot 
component. It seems that the ionization state of Si in the hot compo¬ 
nent is somewhat too high. This may be caused by the fact that this 
spectrum was taken from a large region of the remnant, in which 
in different parts the Si may be in different ionization states. The 
parameters obtained from this spectral fit will be used in the discus¬ 
sion to determine the overall physical parameters of the remnant. 

Our best-fit model differs significantly from results obtained 
previously. The most recent X-ray observations of 3C400.2 were 
performed by |Yoshita et al.||2001[ >, using the GIS instrument on¬ 
board ASCA X-ray telescope. The Nh that we find is quite similar 
to theirs, but they find an acceptable fit for this region using a single 
NEI model with kT = 0.8 keV and log(r) = 10.7 — 11.2 cm -3 
s. We have tried to fit our data with a single NEI model like this, 
but this gives a C-stat / d.o.f. = 606.23 / 229, even when allowing 
multiple abundances to vary. The differences in best-fit parameters 
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Figure 5. A close up of the full spectrum, showing the effects of overioniza¬ 
tion in 3C400.2. The red solid line shows the best fit model without using 
cooling, while the blue line shows the best fit cooling model for this region. 
The bottom panel shows the residuals of the non-cooling model with plotted 
in blue the difference between the cooling and non-cooling model divided 
by the error. This shows that the cooling model fits the spectral features 
around 0.87 keV and 1.4 keV much better than the non-cooling model. 

most likely stem from the fact that the Chandra ACIS instrument 
has a higher spectral resolution than the GIS instrument, allowing 
for better constraints on plasma parameters. 

2.2.1 Overionization 

Although overionization is a common feature of MMRs, it is not 
immediately clear that the plasma is overionized in 3C400.2, due 
to the lack of strong, isolated RRC features. Overionization in this 
remnant has a more subtle presence, which we illustrate in Fig. 15] 
This figure shows a close-up of the full spectrum shown in Fig. [3] 
in the energy range 0.5 - 1.7 keV, where the red line represents the 
best-fit model without overionization, and the dashed blue line the 
best-fit overionized model. The bottom panel of the figure shows 
the residuals of the non-overionized model, with plotted as a blue 
line the two models subtracted and divided by the error. It is clear 
from the figure that the overionized model fits the data much better, 
especially in the 0.7-1.0 and the 1.3-1.7 keV region. In the 0.5-1.7 
keV range, the best-fit single absorbed NEI model has a C-stat / 
d.o.f. = 334/72, a double absorbed NEI model has C-stat / d.o.f. = 
119/64, and the overionized model with the parameters shown in 
table[T|has a C-stat / d.o.f. = 92/70. The overionized model there¬ 
fore fits the data significantly better than any combination of a sin¬ 
gle or double absorbed underionized NEI or CIE model. We have 
performed such a fitting routine where we first attempted a single 
underionized NEI model, then a double underionized NEI model, 
then allowing the initial temperature of the lower kT NEI model 
to vary for every spectrum shown, allowing in every case the abun¬ 
dances to vary only if the fit improved significantly. The Nu was 
always allowed to vary. The C-stat / d.o.f. for the different fitting 
attempts are listed in table [3] 

In the 0.7-1.0keV energy range, the strongest emission feature 
is present at 0.87 keV which is not well-fit by the non-overionized 
model. There are two main spectral emission options which could 
account for the emission feature at 0.87 keV: Fe XVIII emission, 
which has mainly emission lines at 0.77 and 0.87 keV, and the 


Table 1 . The best-fit model parameters of the spectrum covering the whole 
area of the ACIS-I chips. We used the solar abundances from|Anders &| 
|Grevesse] j!989} . Note that the hot ejecta plasma is underionized, while 
the ISM component is oven onized. The errors represent the 1 a confidence 
interval. 


Component 

Parameter 

Unit 

value 

Ejecta 

Ah 

10 21 cm" 2 

e.ost 3 ; 33 


n e n H V 

10 55 cm -3 

o q-i +0.83 
Z - O1 -0.67 


kT e 


q Ofi + 0.30 

O.ob_0.28 


T 

10 10 cm" 3 

9 r\9+0.05 
z - uz -0.06 


Si 


q i 1 +1.23 
0.11 —o.84 


s 


6-09+^ 


Fe 


16-6 til 


Luminosity 

10 31 erg s _1 

68 

ISM 

n e n H V 

10 58 cm -3 

l.lVl'j’g 21 


fcTi 


0 42+ 0,03 
u<4Z -0.02 


fcT e 


0 -m+O-OI 
u.i4_o oi 


T 

10 10 cm -3 s 

26 4“*~ 7,8 


Ne 


Q 40+ 0 04 
u -^ u -0.04 


Luminosity 

10 31 erg s _1 

4.5 


C-stat / d.o.f. 257.84 / 228 


radiative recombination continuum (RRC) of O VIII. Fe-L emis¬ 
sion in general is often not well predicted by plasma models, as 
there are many different emission lines of which the strengths are 
not entirely known or understood (e.g. |Bemitt et al.|2012[ l. Of the 
different ionization states of Fe which produce Fe-L emission, the 
Fe-L lines produced by Fe XVII are often the most prominent in 
the hot plasmas found in SNRs, but they have a weak presence in 
3C400.2. Therefore, for a plasma to show strong Fe XVIII emission 
unaccompanied by Fe XVII emission, the ionization state of the Fe 
plasma needs to be high, with most of the Fe in Fe XVIII-XX ion¬ 
ization states. We have tried to fit this feature solely with a highly 
ionized Fe plasma, but due to other Fe emission lines this always 
resulted in an unacceptable fit, as the higher ionization states of Fe 
produce too much emission around 1.2 keV. The non-overionized 
NEI model attempts to fit the feature at 0.87 keV using a combi¬ 
nation of the Ne IX He — a line and Fe XVIII emission. The Ne 
line has a centroid of 0.92 keV, however, which results in a bad 
fit. The overionized model fits this feature much better due to the 
addition of the O VIII RRC. In the 1.3 - 1.7 keV range, the non- 
overionized model shows large residuals, where it attempts to fit 
continuum emission using lines of Mg XI-XII, which has emission 
lines with centroids of 1.35 and 1.47 keV respectively. The overion¬ 
ized model fits this region much better, using a Ne X RRC at 1.36 
keV, which improves the fit with respect to the non-overionized 
model. 

The above example shows that although the presence of ove¬ 
rionization in 3C400.2 is indeed subtle, the overionized models 
provide a significantly better fit to the data due to the presence of 
RRCs. 


2.3 Spectra 

In this section we apply a cooling model to several regions of the 
remnant, which show significant differences in their plasma proper¬ 
ties. The best-fit models of the different regions are listed in Tab.[3] 
where we also listed the C-stat / d.o.f. of the single NEI and double, 
non-cooling NEI models that we attempted to fit to the data. 
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Table 2. Parameters of the alternative model for region 1. The errors repre¬ 
sent the lcr confidence interval. 


Component 

Parameter 

Unit 

value 

Ejecta 

A h 

10 21 cm -2 

5 61 +0 - 47 

°- Oi -0.45 


n e njjV 

10 54 cm” 3 

2 59+ 028 
z,oy -0.26 


kT e 


q opj+0.89 

o.^U—o 68 


T 

10 10 cm -3 s 

o i 9+0.08 
z - ±z -0.08 


Fe 


015 (fixed) 


Luminosity 

10 31 erg s _1 

4.2 

ISM 

n e n H V 

10 56 cm -3 

7 49+ 213 

‘ —1.41 


fcTi 

keV 

0 77+ 0 - 27 
u -' '-0.15 


kTe 

keV 

0.19+°+ 


T 

10 10 cm” 3 

86.5++° 


Luminosity 

10 31 erg s _1 

1.6 


C-stat / d.o.f. 81/87 


2.3.1 Region 1: southern part of the remnant 

This spectral extraction region coincides with the part of the rem¬ 
nant that is weak in X-ray emission, as is apparent from Fig. [4] 
but is the brightest in terms of optical emission ( |Winkler et al.| 
|1993| . This means that these parts of the remnant either show ra¬ 
diative shocks, where the shock velocity is lower than 200 km s -1 
for which the post-shock plasma cools efficiently, or the plasma 
has otherwise cooled to below 10 G K. This region has two models 
which fit nearly equally well. The first model (C-stat / d.o.f. = 81 / 
87), of which the parameters are shown in Tab. [2] has similar prop¬ 
erties to the model for the full region. It has an ISM component 
with kT e < kTi and an ejecta component with kT e = 3.1 keV. 
Compared to the full model, however, the ejecta component con¬ 
tains no overabundance of Si. The feature in the spectrum around 
1.8 keV is fit by the Mg IX RRC of the cooling component. The 
r = 3.20lg'gg x 10 10 cm”* 1 s combined with the high temper¬ 
ature produces a Si line with a centroid at ~ 1.82 keV, which is 
too high to fit the feature. In addition, the Fe abundance is very 
poorly constrained: when left as a free parameter the best fit value 
raises to ~ 250 times solar abundance, suggesting a plasma con¬ 
sisting purely of Fe. This is not realistic for 3C400.2, especially 
when considering the best fit values of the other regions, and there¬ 
fore we fixed the Fe abundance at 15, close to the value obtained 
for the full region. The spectrum with the best-fit model is shown 
in Fig. [6] 

The second model for this region (C-stat / d.o.f. = 78.2 / 84) 
consists of two cooling NEI components of which the parameters 
are listed in Tab. [3] The coolest component has kT e = 0.06 keV, 
the lowest electron temperature of all regions in the remnant. The 
spectral features around 0.8-0.9 keV (see Fig. [7j are in this model 
fit by a combination of O VII-VIII RRCs at 0.74 and 0.87 keV. 
In addition it shows a Ne IX RRC at 1.2 keV. There is virtually 
no line emission present in this component. The hotter component 
with kT e = 0.53 keV shows the O VIII emission line coupled with 
an RRC, but otherwise again very little line emission. The abun¬ 
dances of both components are mostly solar, with only O being 
overabundant in the hotter component. Contrary to the full spec¬ 
trum and spectra of other extraction regions, there is no super solar 
abundance of Si, S or Fe present in the hot component for this re¬ 
gion. The Fe abundance is sub-solar in the hotter component, while 
in the cooler component the temperatures are such that Fe has an 
ionization state <Fe XVII, which does not show significant emis¬ 
sion in the X-ray band. The fact that there is practically no line 
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Figure 6. The spectrum of region 1, fitted with the model of which the 
parameters are listed in Tab.[2] The red solid line shows the complete model. 
The blue long dashed line shows the low kT NEI component while the 
green short-dashed line shows the high kT component. 

emission from elements with higher mass than Mg is not unex¬ 
pected in a plasma that is cooling rapidly, since the recombination 
and ionization rates of an element depend strongly on its charge. 

The two models are statically almost equally probable. The 
first model fits the 0.8-0.9 keV features slightly better, while the 
second model fits the data better in the 0.4-0.8 keV region. Based 
on the abundances, especially the absence of Si and the high over 
abundance of Fe, we deem the double cooling NEI model slightly 
more realistic. The strong cooling fits well with the spatial coinci¬ 
dence with the optical emitting region in the remnant. 

2.3.2 Region 2: brightest part of the remnant 

This region was extracted from the brightest part of the remnant in 
X-rays (see Fig.[4]l. The spectrum is shown in Fig. [7] top right. The 
best fit model again contains a hot ( kT e = 3.59 keV) underionized 
(t = 1.87 x 10 4 cm -3 s) NEI component with super-solar abun¬ 
dances, plotted as a green dashed line, and a cooler ( kT e = 0.63 
keV) NEI component that is overionized, plotted as a blue dot- 
dashed line. As expected, the best fit model is very similar to the 
model of the full region, since the full spectrum is dominated by 
emission from the brightest part of the remnant. The plasma pa¬ 
rameters of the hot component are identical to the full region within 
the errors, but the Nn is significantly higher at 7.62+°'24 x 10 21 
cm -2 , compared to 6.081q+ x 10 21 cm~ 2 . In addition, the initial 
temperature kT\ = 0.63 keV for the cooler model is somewhat 
higher in this region than for the central region (fcTi=0.42), while 
the electron temperatures are identical within the errors. In general 
the electron temperatures for the cool component are very similar 
throughout the remnant, while the initial temperature shows signif¬ 
icant variation. 

2.3.3 Region 3: The northwestern part of the remnant 

This region was taken from the utmost NW part of the Chandra 
FOV, which overlaps with the radio shell. The spectrum is shown 
in Fig.[7]bottom left. The hydrogen column of the best-fit model of 
this region is 9.7 x 10 21 cm -2 . Although equal within the errors 
to the Nh of region 2, it is significantly higher than the rest of the 
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Table 3. Best-fit model parameters for the spectral regions shown in Fig. [2] The abundances of the unlisted 
elements are fixed at their solar values. The errors represent a ler confidence interval. The C-stat / d.o.f. of the 
single NEI and double, non-cooling NEI models we attempted to fit to the data are listed. 


Region 


Parameter 

Unit 


1 

2 

3 

4 

A h 

10 21 cm' 

-2 

4 76+ 0,45 
<O_ 0.88 

7 co+0.26 

‘ —0.24 

9-72+1 

5 49+ 023 

J -^ y -0.23 

n e n nV 

10 58 cm' 

-3 

0 03 +0 01 
u.uo _ 0 oi 

179++(xl0- 3 ) 

8 - 80 + 7 o( Xl0_3 ) 

0.275+'+ (xlO -3 ) 

fcTi 

keV 


0 53 + 0 - 04 
U -°°-0.04 

- 

- 

- 

kT e 

keV 


0 25 +0 - 04 
u - zo - 0.02 

3-59 + 44 

0 71 +°- 06 

u - / J --0.13 

q 9 q + 0.67 

O.ZO_072 

T 

10 10 cm' 

- 3 s 

1 Q7+2.77 

J--o< _1.87 

9 nq+0.06 

Z.UO-006 

19.58+1 

1-94+iI 

o 



2 67 +0 - 91 
-0.61 

- 

- 

- 

Ne 



< 0.24 

- 

- 

- 

Si 



- 

9 7c:+0-81 

Z -'°-0.57 

1-34++ 

9.53++ 

S 



- 

5-05± 4 ;® 2 

2-98+; 22 

12 . 2+° 9 8 

Fe 



0 62”*” 0,20 
U - OZ -0.24 

12 . 8+1 

0 47+0-28 
-0.47 

31.8+ 8 5 2 

Luminosity 

erg s -1 


2.7 x 10 31 

4.1 x 10 32 

9.53 x 10 31 

1.4 x 10 32 


n e n H U 

10 58 cm -3 

0 16”*~°' 15 
u - ±d -0.06 

q 01 +0.07 
u,o± -0.06 

2.08++ 

0 2 Q +0 07 
u ’ zy -0.05 

fcTi 

keV 

0 21 +0 08 
u - zl - 0.02 

0 63 - * -0 11 
u - oo -0.07 

0.15++ 

r\ q9+0.03 
u ‘ oz -0.04 

kT e 

keV 

0 Ofi”*" 0 01 
u.uo_ 0 .oi 

0-11++1 

O-lO+oi 

0.14++ 

T 

cm -3 s 

9-40+°4q(x 10 10 ) 

5.49++(xlO n ) 

2.18++(xl0 8 ) 

4.63++X10 10 ) 

O 


- 

q n o+0.53 

3 .U 8—0 49 

- 

- 

Ne 


- 

- 

9 H7+0.25 
z,u '-0.51 

0 59 +0 ' 35 
u - oy -0.i5 

Luminosity 

erg s —1 

5.3 x 10 26 

1.6 x 10 31 

7.2 x 10 28 

7.5 x 10 3 ° 

C-stat / d.o.f. 


78.2/84 

118.40/103 

93.68 / 108 

90.57 / 90 

C-stat / d.o.f. single NEI 


140.5 / 87 

226/ 106 

135/110 

253 / 93 

C-stat / d.o.f. double NEI 


115/86 

163 / 102 

120/107 

149/88 


remnant. An increasing Ah towards the NW of the remnant is con¬ 
sistent with the notion of the remnant expanding into a dense ISM 
cloud in the NW, where it is expected that the emission from the 
outer parts of the remnant travel through a larger amount of mate¬ 
rial than the emission from the inner parts. Significant variation in 
foreground absorption cannot be ruled out, however. 

The best-fit model again has a somewhat hotter, underion¬ 
ized NEI component with super-solar abundances combined with a 
cooler, overionized NEI component. However, the electron temper¬ 
ature of the hotter component is significantly smaller at kT e = 0.71 
keV, than those of the full spectrum and region 2 (« 3.75 keV), 
while the r is significantly larger at ~ 2 x 10 11 cm' 3 s, compared 
to t = 2 x 10 4 cm -3 s. And while the abundances are signifi¬ 
cantly super solar for Si and S, the Fe abundance is lower. It seems, 
therefore, that the Fe emission is confined to the central parts of the 
remnant. 

2.3.4 Region 4: The centre of the remnant 

Region 4 was taken from the SW part of the Chandra FOV, which 
is situated roughly in the centre of the total remnant, as can be de¬ 
duced from the ROSAT image. The spectrum of this region has 
again similar parameters as the spectra from the full remnant and 
region 2. It is plotted in Fig.[7]bottom right, where the hot NEI com¬ 
ponent is plotted in a green dashed line, while the blue dot-dashed 
line represents the cooler NEI component. The Ah = 5.49 x 10 21 
cm -2 of the best-fit model for this region confirms the trend of 
increasing Ah towards the NW part of the remnant. The plasma 
parameters of the hot component are equal within the errors to the 


temperatures of the full spectrum and the spectrum of region 2. 
Overall the model for this region confirms the notion of hot ejecta 
being confined to the central, brighter part of the remnant. 


3 DISCUSSION AND CONCLUSION 

From our spectral modelling we obtain the following. The overall 
spectrum is well-fitted by a two component NEI model plasma, of 
which one is tmderionized with a high 0.71 < kT < 3.86 keV 
and super-solar abundances of Si, S and Fe, while the other NEI 
component has a lower 0.06 < kT e < 0.15 keV, is overionized 
and has approximately solar abundances. The central parts of the 
remnant show significantly higher abundances than the outer parts, 
which is apparent in region 1, where no super-solar abundances are 
found. Region 1 coincides with the optically emitting region, and 
shows the lowest electron temperature. Although different parts of 
the remnant show slightly different plasma properties in terms of 
initial temperature kT\ and electron temperature T e , the parame¬ 
ters are consistent with a hot ejecta plasma confined to the central 
part of the remnant, which is surrounded by a rapidly cooling, ove¬ 
rionized swept up ISM plasma. 

3.1 Shocked mass 

Using the parameters listed in Tab.|7] we can estimate of the amount 
of X-ray emitting mass, both in ejecta and ISM, present in the SNR. 
The emission measure of the cool component is 1.2 ± 0.2 x 10 58 
cm' 3 . We use a distance of 2.5 kpc and a spherical volume of the 
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Figure 7. From left to right, top to bottom: spectra of region 1 . region 2, region 3 and region 4. The models are described in the text. In general the green 
dashed line represents the cooler NEI component, while the dash-dotted blue line represent the hotter NEI component. The total model is plotted as a solid red 
line. Note that 


emitting region. The spectral extraction region has a radius of 6.67 
arcmin (4.85 pc), which corresponds to Vtotai = 1-4 x 10 58 d 3 . 5 
cm 3 . If we assume that the hot and the cool component are two sep¬ 
arate plasmas, which both occupy part of the total emitting volume 
and which are in pressure equilibrium, we can get a unique solution 
for the density and shocked mass of the different components. The 
reason is that if Phot = Pcool then also nhotkT ho t = n coo ikT coo i, 
where n is the number density. For n e = 1.2 nH, the number density 
71 h of a component is given by {EM /1.2/V) 1 ^ 2 , so that: 


EM coo i 
1.2 Vtotai (*) 


1/2 


kT c , 


( EM hot \ 
V 1.2Vtotal(l — X) J 


1/2 

kThot, 


( 1 ) 


where EM = n e nnV, Wot = (1 - x)V to tai and H coo i = 
xVtotai- The above equation is equal for x = 0.4, so that 
ricooi = 1.3(d 2 .5) _1/2 cm -3 and rihot = 0.05(d 2 . 5 ) _1/2 cm -3 . 
The respective masses are M coo \ = 8.7(d 2 .s) 5 ' 2 Mg and 

Mhot = 0.46(d 2 .5 ) J/, “ Mg. 

The total emitting volume is uncertain, since the line of sight 
depth of the remnant is unknown and therefore the emitting volume 
might be a factor of two greater. Note that the total mass in the 


remnant is larger than we calculate here, since the Chandra FOV 
covers about half of the total area of 3C400.2, and we only estimate 
the mass of the plasma with T > 10° K. 

The total mass of ejecta and ISM is quite similar to the 
mass found by |Yoshita et al.| l |2001| l, who found a mass of 6.7 ± 
1.2 Mgd^s f° r the whole remnant. This is surprisingly small for a 
mature remnant like 3C400.2, even if the remnant would be located 
at twice the assumed distance. However, it should be noted that this 
is only the X-ray emitting mass, and most of the mass may in fact 
be ‘hiding" in the plasma cooled below 10 6 K. The hot component, 
with super-solar abundances, only makes up a small fraction of the 
mass, and the total ejecta mass seems very low if it is the ejecta 
component of a massive star. However, here it should be noted that 
only the inner regions of a massive star have enhanced metallicity. 
The total metallicity is a function of stellar mass, with stars around 
13 Mg producing only 0.3 Mg of oxygen l |Vink|2012afr . Therefore 
a low mass of ejecta-rich material is consistent with a relatively low 
mass for the exploding massive star. This does, however, also sug¬ 
gest that the cooler plasma, which we designated ISM, may par¬ 
tially or completely consist of the hydrogen-rich envelope of the 
star. An explanation for why, in particular, ejecta material remains 
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hot needs to be addressed by detailed hydrodynamical simulations, 
which should incorporate the effects of the stellar wind of the pro¬ 
genitor. 

The low mass of the metal-rich component is also consistent 
with a Type la origin for the remnant. In that case the cooler compo¬ 
nent could be solely shocked ISM. The enhanced iron abundances 
may indeed hint at a Type la origin, although the total mass of Fe at 
an abundance of 15 times solar is still much lower than the H mass. 
In addition, the association with the HI regions in the NW part of 
the remnant makes a core-collapse origin more likely. To settle on 
the origin of the remnant, it would be important to reconfirm the 
distance estimate of 2.5 kpc by jGiacani et al. ( } 1998^ , as a larger 
distance estimate would favour a core collapse origin. It would also 
be helpful to identify a stellar remnant in 3C400.2. The Chandra 
image does not show any evidence for a bright point source that 
could be the cooling neutron star. 

3.2 Evolutionary scenario 

The above densities suggest a hot, metal enriched, tenuous plasma 
surrounded by a dense, cooler plasma which is cooling rapidly. This 
is consistent with the shell-like density structure expected from a 
Sedov evolutionary scenario. This is not the only MMR remnant 
in which such a temperature and density gradient is observed as 
|Kawasaki et al.|{2002} also find a two temperature best-fit model 
and overionization in IC 443. They find that the central ejecta re¬ 
gion has a temperature of 1 keV, compared to 0.2 keV in the outer 
layers. However, overionization is present only in the central ejecta 
rather than in the outer layers, while we find overionization in the 
outer layers and not in the ejecta. 

In the introduction we mentioned three different evolutionary 
scenarios that might explain the centrally peaked X-ray emission: 
the evaporating cloudlet scenario (White & Long 1991), the |Cox| 
et al. (19991 scenario with high surrounding density including ther¬ 
mal conduction, and the relic X-ray emission scenario. Our obser¬ 
vations show that the plasma is best explained by a low density, hot 
interior surrounded by a high density, lower temperature plasma. 
The total X-ray emitting mass is relatively low, for a mature rem¬ 
nant. But this is likely an indication that most of the remnant mass is 
not emitting in X-rays. In addition, overionization is only present in 
the cool plasma, suggesting that it is cooling more efficiently than 
the hot plasma. Our results are perhaps most naturally explained 
by the simplest evolutionary scenario of the three: the relic X-ray 
emission scenario ( |Rho & Petre|1998j >. As a result of the high sur¬ 
rounding ISM density, the outer layers have cooled below temper¬ 
atures capable of emitting in X-rays and the interior is still hot but 
has a low density and is therefore not cooling efficiently. We can¬ 
not rule out the presence of thermal conduction, but we do not find 
evidence for overionization as a result of rapid cooling for the hot 
centrally confined plasma. Moreover, the X-ray emitting plasma in¬ 
side the remnant is clearly not isothermal, as indicated by the model 
of |Cox et al.| l |1999) . As a final note, thermal conduction has mainly 
been introduced to models explaining the evolutionary scenarios of 
MMRs based on the then current observations of generally lower 
spectral and spatial resolution, which showed little to no tempera¬ 
ture gradient in the remnants. However, more modeling is needed 
to understand whether local density alone determines if a remnant 
will evolve into a MMR. or whether some other conditions, such as 
pre-supernova evolution or ejecta structure, are important as well. 

The overionization of thermal plasmas can quite naturally oc¬ 
cur in MMRs. The high initial ISM density allows the plasma to 
reach CIE on a timescale smaller than the age of the remnant, after 


which a combination of adiabatic and radiative cooling can make 
the cooling rate of the plasma higher than the recombination rate, as 
shown in Broersen e t al.| (20ll). The higher surrounding ISM den¬ 
sity of MMRs might then be the determining factor for the occur¬ 
rence of overionization compared to non-MMRs, as already noted 
in |Vink| ( j2012b) . Indeed, all remnants cool adiabatically and by ra¬ 
diation, but not all remnants expand in a high enough ISM density 
to reach CIE and then overshoot to overionization. 


4 SUMMARY 

We have analysed an archival Chandra observation of the mixed 
morphology remnant 3C400.2. Our results can be summarised as 
follows: 

• The plasma of the mixed-morphology SNR 3C400.2 is best 
fitted by a combination of a hot, underionized plasma with low 
density, and a cooler, overionized plasma with high density. To our 
knowledge, this is the first evidence for a combination of an ove¬ 
rionized outer shell surrounding an ejecta-rich, underionized inner 
region in an MMR. 

• The hot plasma shows significant overabundances of Fe, Si 
and S, suggesting an ejecta origin, with Fe enhanced in the central 
part. 

• Overionization is significantly present in all parts of the rem¬ 
nant covered by the Chandra field of view. 

• The X-ray emitting masses of the plasma components are 
Mcooi = 8.7(d 2 .s) 5/2 M 0 and M hot = 0.46(d 2 . 5 ) 5/2 M 0 . 

• This low overall mass suggests that most of the X-ray emitting 
mass is from mix of metal-rich and hydrogen-rich (envelope) ejecta 
from a not too massive core collapse supernova, or the remnant has 
a Type la origin. 

• The observations are best explained by a scenario in which the 
centrally peaked X-ray emission is caused by a hot, metal enriched, 
tenuous plasma. Due to the high surrounding ISM density the outer 
parts of the remnant have cooled efficiently towards a temperature 
below which they do not radiate in observable X-ray emission. 

• The overionization can be naturally explained by efficient 
cooling due to a high ISM density in combination with adiabatic 
expansion. 
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